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Abstract 

The differential systems satisfied by orthogonal polynomials with arbitrary semiclassical measures sup- 
ported on contours in the complex plane are derived, as well as the compatible systems of deformation 
equations obtained from varying such measures. These are shown to preserve the generalized monodromy 
of the associated rank-2 rational covariant derivative operators. The corresponding matrix models, con- 
sisting of unitarily diagonalizable matrices with spectra supported on these contours are analyzed, and it 
is shown that all coefficients of the associated spectral curves are given by logarithmic derivatives of the 
partition function or, more generally, the gap probablities. The associated isomonodromic tau functions 
are shown to coincide, within an explicitly computed factor, with these partition functions. 
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1 Introduction 

The partition function for Hermitian random matrix models with measures that are exponentials of a 
polynomial potential was shown in ^ to be equal, within a multiplicative factor independent of the 
deformation parameters, to the Jimbo-Miwa-Ueno isomonodromic tau function |11| for the rank 2 linear 
differential system satisfied by the corresponding set of orthogonal polynomials. The results of @] were 
in fact more general, in that polynomials orthogonal with respect to complex measures supported along 
certain contours in the complex plane were considered. These may be viewed as corresponding to unitarily 
diagonalizable matrix models in which the spectrum is constrained to lie on these contours. 

The purpose of the present work is to extend these considerations to the more general setting of com- 
plex measures whose logarithmic derivatives are arbitrary rational functions, the associated semiclassical 
orthogonal polynomials and generalized matrix models. By also including contours with endpoints, the 
latter viewed as further deformation parameters, the gap probability densities are included as special 
cases of partition functions. 

To place the results in context, we first briefly recall the main points of @], restricting to the more 
standard case of Hermitian matrices and real measures. Consider orthogonal polynomials 7r„(x) € 
L^(M, e"'' ^^^^^(ix) supported on the real line, with the measure defined by exponentiating a real poly- 
nomial potential 

V{x) = (1-1) 
j=i 

(Here we assume V{x) is of even degree and with positive leading coefficient, although these restrictions 
are unnecessary in the more general setting of 4 .) The small parameter ft is usually taken as 0{N^^) 
when considering the limit N oo. 

Any two consecutive polynomials satisfy a first order system of ODE's 

dx \ TTn{x) J ^ V ^n{x) J ' ^ ^ 

where V^ix) is a 2 x 2 matrix with polynomial coefficients of degree at most d — 1 — Aeg{V' {x)). 
The infinitesimal deformations corresponding to changes in the coefficients {tj} result in a sequence of 
Frobenius compatible, overdetermined systems of PDE's 

f^ir{"i-tA-T^A-){^rtA j-i,...,d, (1-3) 



dtj \ TTn{x) J ri, \ J y 7r„(a;) 
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where the matrices Tnj{x) are polynomials in x of degree J which satisfy the compatibility conditions 



d d 

fi-TH Tn,j{x),h- Vnix) 

at,] ax 



(1-4) 



It follows that the generalized monodromy data of the sequence of rational covariant derivative operators 
^(E ~ T^nix) are invariant under these deformations, and independent of the integer n. 

This is a particular case of the general problem of rational isomonodromic deformation systems I I . . 
An important role is played in this theory by the isomonodromic tau function associated with any 
solution of an isomonodromic deformation system. This function on the space of deformation parameters 
is obtained by integrating a closed differential whose coefhcicnts are given by residues involving the 
fundamental solutions of the system. The main results of @| were the following. First, the coefficients of 
the associated spectral curve, given by the characteristic equation of the matrix !)„ (x) , can be obtained 
by applying certain first order differential operators with respect to the deformation parameters (Virasoro 
generators) to hi{Zn{V)), where the partition function Z„(V) of the associated n x n matrix model is 



Zn{V) 



dM exp 



--Try(A/)j 



(1-5) 



Second, this partition function is equal to the isomonodromic tau function up to a multiplicative factor 
that does not depend on the deformation parameters 



r,f^ = Z„(y).F„ 



(1-6) 



The present work generalizes these results to the case of measures whose logarithmic derivatives are 
arbitrary rational functions, including those supported on curve segments in which the cndpoints may 
play the role of further deformation parameters. The latter are of importance in the calculation of 
gap probabilities in matrix models |17j since these may, in this way, using measures supported on such 
segments be put on the same footing as partition functions 6 . A Frobenius compatible system of first 
order differential and deformation equations satisfied by the corresponding orthogonal polynomials is 
derived (Propositions 13.21 13. 4|) and the coefficients of the associated spectral curve are again shown to 
be obtained by applying suitable Virasoro generators to ln(Zjv(T^)) fTheorems 14. llH?^ . A formula that 
generalizes (|l-6|l is also derived fTheorem lS.in : 



JM 



(1-7) 



where the factor .7-"„(y) is an explicitly computed function of the deformation parameters determining 
y, which can in fact be eliminated by a making a suitable scalar gauge transformation. 

The results of Theorems 14.11 14.21 give a precise meaning, for finite n, to formulae that are usually 
derived in the asymptotic limit n — > oo (?in ^ 0(1)) through saddle point computations, relating the free 
energy to the asymptotic spectral curve. It is well known ^ that the free energy in the large n limit is 
given by solving a minimization problem (in the Hermitian matrix model) 



^0 



lim h lnZ„ = min 



V{x)p{x)dx 



p{x)p{x') In |a; — x'l 



(1-8) 



giving the equilibrium density pcq for the eigenvalue distribution. If, for instance, the potential is a real 
polynomial bounded from below, it is known jHj that the support of the equilibrium density is a union of 
finite segments / C M. The density poq is obtained from the variational equation 



2V 



Pcq{x)dx 



V'{x) 



(1-9) 
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(within /), and is related to the resolvent by 









-I. 







dx^^ , zeC\l . (1-10) 

X — z 

In terms of this, the spectral density may be recovered as the jump-discontinuity of uj(z) across J, and 
all its moments are given by 

1 f x'^ z"^ 

dtjJ^o = - lim h(TrM-^) = / dx — p^Jx) = - res — u;{z)dz . (1-11) 

J n^oo J I J z=oc J 

The function y = —uj{x) satisfies an algebraic relation given by 

=yV'{x)+R{x) , (1-12) 
where R{x) is a polynomial of degree less than V'{x) that is uniquely determined by the consistency of 

(inni and (nnai- 

The point to be stressed here is that this asymptotic spectral curve should be compared with the 
spectral curve of Theorem 14.11 given by the characteristic equation H4-29|l . which also contains all the 
relevant information about the finite n case. In the n —^ oo limit, logarithmic derivatives of the partition 
function are expressed in l|l-ll() as residues of the meromorphic differentials z'^ydz on the curve. The 
same formulae are shown in Theorem 14.21 to hold as exact relations for the finite n case if we replace 
the "asymptotic" spectral curve by the spectral curve given by the characteristic equation of the matrix 
■Dnix). 

The paper is organized as follows. In section |21 the problem is defined in terms of polynomials or- 
thogonal with respect to an arbitrary semiclassical measure supported on complex contours, and the 
corresponding generalized matrix model partition functions. In section [2. 21 the recursion relations, differ- 
ential systems and deformation equations which these satisfy are expressed in terms of the semi-infinite 
"wave vector" formed from the orthogonal polynomials. In sectionOthe notion of "folding" is introduced 
and used (Propositions 13 . 4|1 to express the preceding equations as an infinite sequence of compatible 
overdetermined 2x2 systems of linear differential equations and recursion relations satisfied by pairs of 
consecutive orthogonal polynomials. In section 0] the results of folding arc used to express the spectral 
curve in terms of logarithmic derivatives of the partition function and it is shown that the n — > oo relation 
between the free energy and the spectral curve is also valid as an exact result for finite n. In sectional 
the definition of the isomonodromic tau function 11 is recalled and it is computed by relating it to the 
spectral invariants of the rational matrix generalizing P„ (x) in Hl-2|) . These invariants are shown to give 
the logarithmic derivatives of the tau functions in terms of residues of meromorphic differentials on the 
spectral curves through formaulae that are nearly identical to those for the partition function, This leads 
to the main result, Theorem 15. II which gives the explicit relation between Z„ and r^*^. 



2 Generalized orthogonal polynomials and partition functions 

2.1 Orthogonality measures and integration contours 

Given a measure on the real line, the associated orthogonal polynomials are those that diagonalize 
the quadratic form associated to the corresponding (complex) moment functional; i.e., the linear form 
obtained by integration with respect to the measure. 

C:C[x] ^ C. 

p{x) i-^ C{p{x)) = J^p{x)diJ.{x) . 



4 



A natural generalization consists of including moment functionals that are expressed by integration along 
more general contours in the complex x plan, with respect to a complex measure defined by locally 
analytic weight functions that may have isolated essential singular points and complex power-like branch 
points. 

We thus consider linear forms on polynomials given by integrals of the form 



C{p{x)) — J p{x)ii{x)dx 

K 

Vix):=J2Tr{^) ^ (2-2) 



r=0 



where 



to, J 



To{x) := to,o + ^ —^x" 



dr 



tr,J 



Tr{x) := ''" rj - Ufl \tl{x - Cr) 

j^^ J{x-Cry 

do-1 K dr + l 

- hd,, \r,^x{x) = V'{x) = to,jx'-' ( ' (2-3) 

J=l r=l J=l ^'■> 

and the symbol denotes integration over linear combinations of contours on which the integrals are 
convergent, as explained below. This class of linear functionals is sometimes referred to as semiclassical 
moment functionals pi ll4[[T5) . We consider the corresponding monic generalized orthogonal polynomials 
Pn{x), which satisfy 

P„{x)pm{x)n{x)dx = KSnm ■ (2-4) 



If all the contours are contained in the real axis and the weight is real and positive, we reduce to the 
usual notion of semiclassical orthogonal polynomials. The small parameter h introduced in (|2-2|l is not of 
essential importance here; it is only retained in the formulae below to recall that, when taking the large 
n limit, it plays the role of small parameter for which hn remains finite as n — > oo. 

To describe the contours of integration, we first define sectors Sr'^\ r — 0, . . . K , k — 1, . . . d^.. around 
the points for which > 1 (cq := oo) in such a way that 

5R(F(a:))^^_^ +00 . (2-5) 

The number of sectors for each pole in V is equal to the degree of that pole; that is, do for the pole at 
infinity and dj. for the pole at c^. Explicitly 

^r-j.eC; ~ - ^ < arg(x) < " ^'^^J^^"-'^" ^ + ^ } , (2-6) 

[xe C; ' '^^ + ^^^i^r,d,.)~- ^ ^^^(^ _ ^ 2fc. + argft.,.J + || 

fc = 0, r^l,...,K. 

These sectors are defined in such a way that approaching any of the essential singularities of /i(x) (i.e. a 
Cr such that dr > 0) within them, the function fi{x) tends to zero faster than any power. 
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2.1.1 Definition of the boundary- free contours 

The definition of the contours follows ^Hl (see fig. 1). 

1. For any Cr for which there is no essential singularity in the measure (i.e. dr = 0), there are two 



(a) For the c^'s that are branch points or poles in fj, (i.e., tr,o ^ N), we take a loop starting at 
infinity in some fixed sector encircling the singularity and going back to infinity in the 
same sector. (Note that if Cr is just a pole; i.e., tr,o G — N+, the contour could equivalently be 
taken as a circle around c^.) 

(b) For the Cr's that are regular points (t^^o G N ), we take a line joining Cr to infinity, approaching 
oo in a sector s'f'^ as before. 

2. For any Cr for which there is an essential singularity in fi (i.e. , cZi- > 0) we define dr contours 

(r^ . . . (r) A . ± 

starting from Cr in the sector 5*^ and returning to it in the next sector Sj^_^-^. Also, if trfi ^ Z, we 
join the singularity Cr to cxd by a path approaching oo within one fixed sector 5'^'^''. 

3. For Co '■— CXD, we take do ~ 1 contours starting at cq in the sector and returning at co in the 
next sector S^^^. 

Note that, with these definitions, the integrals involved are convergent and we can perform integration 
by parts. Moreover, any contour in the complex plane for which the integral of ii{x)p{x)dx is convergent 
for all polynomials p{x) is equivalent to a linear combination of the contours defined above, no two of 
which are, in this sense, equivalent. 

2.1.2 Definition of the hard-edge contours 

We also include some additional contpours in the complex plane {m)j}j^i^,,,^L, starting at some points 
Gj, j — 1 . . . L and going to cxd within one of the sectors S^^"^ . These could be viewed as corresponding to 
additional points in 1(b) for which both dr = and ir,o — 0, but we prefer to deal with them separately 
since integration by parts on these contours does give a contribution. 

In total there are S := do + ^r=ii'^r + 1) boundary-free contours ag, ,£ = 1, . . . , S and L hard-edge 
contours m^, /i = 1, . . . , L. The moment functional is an arbitrary linear combination of integrals taken 
along these contours 



Note that, by taking appropriate linear combinations of the contours, we could alternatively have had 
contours consisting of finite segments joining the points aj. 



subcases: 




(2-8) 



6 



Figure 1: The types of contours considered in the x Riemann sphere P^. Here 
we have ci with di = 3 and C2 with d2 = 0, t2,o ^ Z (logarithmic singularity 
in the potential), C3 with da — 0, ^3,0 G N and the degree of the potential at 
infinity co = 00 is do = 5. The essential singularity in fi at ci is of the form 
exp {x — ci)~^ and there is also a cut extending from ci to 00 if ti^ ^ Z. The 
point C2 is a branch point of fi{x) since ^2,0 ^ and the cut extends to infinity 
"inside" the contour (as shown here). If it were a pole (^2,0 G — N+), the contour 
would be replaced by a circle around it. The point C3 is a regular point with 
^3,0 € , and the contour extending from it to infinity is no different from the 
ones starting at the regular points 01,02,03. The latter are the "hard-edge" 
segments joining the points oi , 02 and 03 to 00 within one of the sectors S^'^ . 



2.2 Recursion relations, derivatives and deformations equations 

2.2.1 Existence of orthogonal polynomials and relation to random matrices 

Recall j2j that orthogonal polynomials satisfying H2-4|l exist provided all the Hankel determinants formed 
from the moments are nonzero: 



Since the A„(>f)'s are homogeneous polynomials in the coefficients Xj, the zero locus excluded by (|2-9(l 
is of zero measure (in the space of Xj^s), and hence "generically" the conditions (|2-9|l are fulfilled. The 
development to follow will in fact only involve orthogonal polynomials up to some arbitrarily large fixed 
degree, say N, and hence the conditions A„(>tf) =0 ,n < N — 1 determine a Zariski closed set in {xj}, 
(and a closed set of measure zero in the space of coefficients of V) . 

The orthogonal polynomials considered here are related to models of unitarily diagonalizable random 
matrices M e gl{n,C) with spectra supported on the contours defined above. More specifically we have 
the partition function 




(2-9) 
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where 



: J X 

= n\An{K, V)^nlY[ hj , (2-10) 

3=0 

Cn := 7 ^- ^ (2-11) 

is the inverse of the U{n) group volume, and 

A(ai) := Yl(x, - xj) (2-12) 

is the usual Vandermonde determinant. The notation spec{M) £ x in the first integral just means 

M = UDU^, D :=diag(a;i,...a;„) , UeU{n), (2-13) 

where the eigenvalues {2:1, . . of M are constrained to lie on the contours entering in J^. 

In particular, as in the standard case, the orthogonal polynomials may be shown to be equal to the 
expectation values of the characteristic polynomials in such models 

Pn{x) = (dct(a:I-M)) 

= ± f dxi - - [ dx„ f[ix - x,)Aixfe-i^"-^ , (2-14) 

and all correlation functions between the eigenvalues may be expressed as determinants in terms of the 
standard Christoffel-Darboux kernel formed from them 

n— 1 

i^„(x,y) /^P.WP.(2/)e"*^''^"^+''^'" • (2-15) 

More precisely, this is valid when there are no "hard-edge" contours present. Inclusion of the latter 
however allows one to interpret these determinants as certain conditional correlators, known as "Janossy 
distribution" correlators [H], giving the probability densities for a certain number of eigenvalues to lie at 
given locations within the complementary part of the support, while the remaining ones lie within it. 
The partition function Z„ in this case can be reinterpreted as the corresponding gap probability [HI E| ■ 

2.2.2 Wave vector equations 

We now define the normalized orthogonal polynomials 

■Kn{x) ■■= -^Pn{x) (2-16) 

and what will be referred to as the " orthonormal quasi-polynomials" 

i^„(x) :=7r„(x)e-A^(^), (2-17) 

satisfying 

'lpn{x)i^m{x)dx 5mn- (2-18) 
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From the former, we form the semi-infinite "wave vectors" 

n(a;) := (x), . . . , ^„(a;), . . .]* . (2-19) 
As in the theory of ordinary orthogonal polynomials, we have 

a;n(a;) = QIl{x) , (2-20) 
where Q is a symmetric tridiagonal semi-infinite matrix with components 

Qij = 7j'^ij-i + Mi] + Ji6ij+i, ij eN , (2-21) 
defining a three term recursion relation of the form 

XTTj (x) = -fj + lTTj + l (x) + PjTTj (x) + 7j7rj_l (x) . (2-22) 

Now introduce semi-infinite matrices P, Ai,Cr, T^^j such that 

hdjl{x) = Pn(x) (2-23) 

hdaTl(x) = AJl{x) , i = 1, . . . , L (2-24) 

hdcUix) = an(x) , r = l,...,K (2-25) 

hdt,jU{x) = Tr,jll{x) , r^Q,...,K,J = Q,...,dr. (2-26) 

Their matrix elements are determined simply by integration 

Xnrn= / {hdlT n{x)) T^mix) ^.{x)dx , (2-27) 



where d denotes any of the derivatives d^^ dai, dc^, dt^ , above for which X becomes the corresponding 
matrices P, A,, Cr or T^,,/ on the RHS of - 

Remark 2.1 Such wave vectors and associated deformation equations have been studied in many previous 
works relating orthogonal polynomials, matrix models and integrable systems (see, e.g. |2], |18|'). However, 
considerations of the deformation theory have mainly been within the formal setting, with the potential V{x) 
replaced by some initial value, Vo{x), plus a perturbation consisting of an infinite power series with arbitrary 
coefficients, without regard to domains of convergence. Results obtained in this formal setting cannot be directly 
applied to the study of isomonodromic deformations, where the local analytic structure in the neighborhood of a 
number of isolated singular points is of primary interest. 

For any such semi-infinite square matrix X, let Xq, X^, X- denote the diagonal, upper and lower 
triangular parts, respectively, and let 

X_o := ^^0 + . (2-28) 

Proposition 2.1 The matrices P, Ai, Cr and T^j are all lower semi-triangular (with P strictly lower 
triangular) , and are given by 



P = ^'(Q)-o - E = ^'(Q)- - E(^')- (2-29) 

i=l i=X 

A, = ;l>^,(^(a,)^*(a,))-o (2-30) 

Cr = -J2trj{Q-Cr)Zt\ r = \,...,K (2-31) 



J=0 
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To.o = ^1 (2-32) 
To,j^jQio, J=l,...,do (2-33) 

Trj = j{Q-Cr)Zi, r = l,...,K, J^l,...,dr (2-34) 
T,.o = -ln(0-c,)_o, r^l,...,K. (2-35) 
where [Q — Cr)~^ and \\i{Q — Cr) are defined by the formulcE 

(Q _ ,^y^J^ I ^pM^^(,)d, (2-36) 

\tl{Q - Cr)nm '-^ / ln(z - )7r„ (a;)7r„ (z)/i(z)d2; . (2-37) 

T/ie diagonal matrix dements for each of the above is given by the formula 

X,, = -^d{\nh,) , (2-38) 

where d — dx, da^, dc^ and dt^ , , respectively, for X — P, Ai, Cr and T^^j. In particular, they vanish for 
P , which is strictly lower triangular, and hence 

L 

V'iQ)jj = hJ2 ^^^3{a^f ■ (2-39) 

i=l 

Proof. We make use of the orthogonality relations 

U{x)U\x)^i{x)dx = 1 . (2-40) 

Eqs. H2-31|l - l)2-35|l are obtained as follows. Consider a deformation d with respect to any of the above 
Cr's or ir.j's and denote by X the corresponding matrix; then 

na^„(x) = hd^^^ {hd\n{h^)) 7r„(a:) + -^hdpn{x) (2-41) 

V i^n ^ V tin 

= — - {hd\Ti(hn)) 7r„(a;) + lower degree polynomials , (2-42) 

since the polynomials 7r„ are monic. It follows that the deformation matrix X is lower semi-triangular. 
On the other hand, differentiating eq. (|2-40|) gives 

= / (;iann*+nnan*)/i(a;)dx+ / ^^%dy.(x)dx = X + X* + / YlYl^hdy.(x)dx (2-43) 
J ^ J ^ J ^ 

Applying the operators for each case to ji as defined in (|2-2|l and using eq. H2-20|) then gives the result. 

Now consider the deformations of the endpoints ai of the "hard-edge" contours. Differentiating l|2-40|l 
gives 

= hda, j UH*^l{x)dx = -}^^hIl{a^)Il\a^)^l{a^) + j [{hda,Il)U* + UfiBaJl*] fi{x)dx 
= -hx,U{a,)U\a,)tiia,) + A, + A\, (2-44) 
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where 

Ai:=h I daJlU^ii{x)dx , (2-45) 

It follows that 

A, = n>^,(nn*)_oL^„^ , (2-46) 

proving eq. H2-30|l . and also that 

{A{)r,n = -^da, HK) = ^^^^l{a^) . (2-47) 

To determine the matrix P, note that it is strictly lower triangular and 

-h^{A, + Af) ^ h f (n' n* + n n'* + n n*a^ hi ^(x)) ^{x)dx 
1=1 "'^ 

= / {hn'n* + hnu'* -v'{x)nn^) ^{x)dx = p + -v'iQ). (2-48) 

This imphes that 

L L 

= V^'(Q)-o - E ^» = ^'(Q)- - E(^')- ■ (2-49) 

1=1 i=l 

This last equality follows from H2-39I) , which, in turn, follows from integration by parts in the definition of 
y (ff)nn- It may be seen as a consequence of the invariance of the partition function under an infinitesimal 
change in the integration variables Xj — *■ a;^ + e in I|2-1U|1 : i.e., translational invariance. 

From H2-38|l and (|2-1U|I follows a relation between the diagonal elements of the deformation matrices 
and the logarithmic derivatives of the partition function that will be very important in what follows. 
Define the truncated trace of a semi-infinite matrix X to be 

ji-i 

Tr„X:=E^.. ■ (2-50) 

Corollary 2.1 For d = da^, dc^, and dt^ j, 

hdlnZn^ -2TTnX , (2-51) 
with X ~ Aj, Cr and Tr.j, respectively. For the cases dc^ and dt^j, 

ft91nZ„ ^Tr„9y(g) , (2-52) 

while for the da^ 's we have 

L 

J2 l^^^a, ln2„ = -hV'{Q)nn (2-53) 

Proof. The first of these relations follows from H2-38() and (|2-10|l directly, the second from the explicit 
expressions for the deformation matrices (|2-31|l - l|2-35|l and of the potential V{x), and the third is a 
restatement of the H2-39II (translational invariance). 
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Corollary 2.2 The compatibilty conditions 

[G,n]=0, (2-54) 
are satisfied, where G, TL are any of the following operators 

hda,-A„ hdt^,-Tr,j, hdc^-Cr, hd^ - P, x-Q (2-55) 
and r = 0, . . . K , J = 0, . . . dr- 

Proof. This follows immediately from the fact that the orthogonal polynomials entering in eqs. H2-23|l - 
(|2-26|l are linearly independent. 

Remark 2.2 Note that 

[hd^ -P,x~Q]=0 (2-56) 
is just the string equation, while the other compatibility conditions involving x — Q imply the Lax equations: 

hda,Q = [A,,Q], hdt^ jQ = [Tr,J,Q], hdc^Q = [Cr,Q] , (2-57) 

showing that the spectrum of the matrix Q is invariant under these deformations. 

2.2.3 Wave vector of the second kind 

We now consider solutions of the second kind, 



(l)n{x) -.^ e'^^^''' / ^ 



-dz. , (2-58) 

X ~ z 

which may be combined to form the components of a wave vector of the second kind 

*(x) := [Mx),M^),---,M=c),---Y ■ (2-59) 

Denote by 

X 

the semi-infinite square matrix with elements 

■V'{x)~V'iQ)\ f .^iv-f,, ,^,V'{x)-V'{z) 



x-Q 



) = / dze-*^(-)^„(z)7r,„(z)^^^^M_ZM , (2_61) 

/ nm J >^ X Z 



and define U{x) to be the semi-infinite column vector (with only its zcroth component nonvanishing) 
given by 

{U{x))n v^e*^(-)<5„,o. (2-62) 

The following lemma gives the effect of multiplication of 11 (x) by x and of application of hdx to it. 
It may be deduced immediately from eqns. I|2-2()|l and (j2-23|l . apphed inside the integral, together with 
integration by parts. 

Lemma 2.1 

x*(x) = Q*(x) + U{x) (2-63) 

L l-{Vix)-V{a^)) 

hd^^ix) = P^ix) + VqV'{x)U{x) + hY^ X, n(a,) , (2-64) 
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The next proposition, which is similarly verified, gives the effects of the above deformations on the 
wave vector of the second kind. 



Proposition 2.2 

hdt^^^{x 



gJ-(V(x)-y(aO) 



r,J- 



.1=0 



X - flj 

(g - c^-j-i - jx - cQ- j-i 

Q-x 



-Uix), r = l,...,K, 



To,j*(a;) 



Tr,.Mx) 



Q-x 



-U{x), quadJ = 1, . . . , do, 



(2-65) 
(2-66) 
(2-67) 



{Q -~ Cr) - {x- Cr) 
Q ~ X 



Tr,o*(a;) H — U[x), r = l,...,K. 



U{x) , J=l,...,dr, r^l,...,K, (2-68) 

(2-69) 



The content of eqs. (|2-25|l - (|2-26ll and (|2-66|I - H2-69|I may be summarized uniformly as follows. Let 
v{x) be any function that is analytic at each point of the contours except, possibly, the points c^, and 
for which the following integrals are convergent: 



v{Q)r 
{'VQv{x))r 



u(z)7r„(z)7rm(z)e R^f^'^dz = / v{z)^niz)^miz)dz 



v{x) - v{Q) 
x~ Q 



v{x) — v{z) 



1pniz)i>miz)dz. 



(2-70) 



Define the deformation matrix under the infinitesimal variation of the potential V{x) i-^ V(x) + v{x) to 
be 

X, := v{Q)_o . (2-71) 

Then the two infinite systems 



dyU{x) := X^U{x) 

Sy^{x) := Xy^{x) + Wqv{x)U{x) 



(2-72) 
(2-73) 



describe the infinitesimal deformation of the orthogonal polynomials and the second-kind solutions under 
such infinitesimal variations of the potential. 
Equivalently, define the 2 x oo matrix 



T{x) ■= [n(a;),*(a;)] 



(2-74) 



In terms of T{x), all the recursion, differential and deformation equations (|2-20|l . (I2-23|I - (|2-26II and 
(|2-63|I - H2-69|I may be expressed simultaneously as 



X — ai 



xT = Qr+ (0, U) 

hd^T = pr + ( 0, VqV'u + hY^ 

\ i=l 

5^T = X^T + {Q,WqvU) 

hdaT = A,r - 0, hx, n(a,) 

\ X- ai 



U{a,) 



(2-75) 
(2-76) 
(2-77) 
(2-78) 



where v signifies any of the infinitcsilmal deformations of the potential hda, hdt^ ,V {i — 1, . . . L, r 

0,...K, J=l,...,dr). 
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3 Folding 

3.1 n-windows and Christoffel-Darboux formula 

Let in be the oo x 2 matrix that represents the injection of the 2-dimcnsional subspace spanned by the 
(rt — l,n) basis elements into the (semi-)infinite space corresponding to the components of 5* or Its 
matrix elements are thus: 



{in)]k = 5k,l5j^n-l + 4,2(5j,n, j = 0, 1, 2, . 



fc = 1,2. 



(3-1) 



Let denote its transpose, which is the corresponding projection operator. The n-th 2x2 block (or 
"window") of F is then given by: 



Tnix) := ilr = 



(3-2) 



By "folding" the infinite recursion and differential-deformation equations H2-20|l . (|2-23() - (|2-26|l . H2-63|l - 
(|2-69ll . (|2-75|I - H2-78|I . we mean the corresponding sequence of recursion relations, ODEs and PDEs satis- 
fied by the r„(a;)'s. To derive these, a form of the Christoffel-Darboux identity for orthogonal polynomials 
will repeatedly be used. Let IlJ^ denote the semi-infinite square matrix whose only nonvanishing entries 
are I's on the diagonal in positions to n (i.e. the projection onto the first n+1 components) 



(n°). 



Let 



6ij if < i, j < n 
otherwise .quad 

-1 

1 



be the standard 2x2 symplectic matrix, and let 

denote its projection onto the 2x2 subspace in position (n — 1, n). 

Proposition 3.1 The following extended Christoffel-Darboux formulfB are satisfied: 
{x-x')T''{x)Ul^,T{x') 

^ ( \ 

= 7nr„ {x)aTn{x') + y ^^^(^) g 1 ^_ _ j 

-7„r {x)T,nT{x')+ \ ly^^^ i(^v(x)+vix')) f ^-~hv{z)(_i h.-\dz 

\ J yc \x — z x'—z j i 

Equivalently, in components, 



{x - x') ^ TTj{x)TT-i{x') = 7„ [7r„(x)7r„_i(x') - 7r„_i(x)7r„(a;')] 

n-l 

(x - a;') ^ Trj{x)(l)j{x') = 7„ [7r„(x)0„_i (x') - 7r„_i(x)(/)„(a;')] 



.iV(x') 



(3-3) 
(3-4) 
(3-5) 

(3-6) 
(3-7) 

(3-8) 

(3-9) 
(3-10) 



n-l 



(X - x') ^ (f>j{x)(l)jix') = 7„ [(f)n{x)(f)n-l{x') - (/)„_i (x)(/>„ (x')] 
j=0 



14 



^l{Vix) + V{x')) 



and evaluating at x — x' gives 



detr„(x) = TTn^i{x)(pn{x) ~ (/)„- 1 (x)7r„ (x) = LeR^(^) 

In 



(3-11) 



(3-12) 



Proof. Eq. 1)3-91) is the standard Christoffel-Darboux relation for orthogonal polynomials. The extended 
system may be derived as follows. Multiplying the expression T{x)'n'^_iT{x') by {x — x'), and applying 
the relation (|2-75|) with respect to both x and x' gives 



dz 



(3-13) 
(3-14) 



(x - x')T^{x)K_,rix') = r^(x) (gn°_i - n°„_,Q) r(x') 

/ -ei^i-') 

The result 1)3-8(1 is obtained by substituting the following identity, which holds for any tridiagonal sym- 
metric matrix of the form Q 

QUn-i - n"_iQ = -fnincrin ■ (3-15) 



3.2 Folded version of the deformation equations for changes in the potential 

Under infinitesimal changes of the parameters in the potential V and the end-points of the "hard-edge" 
contours, the wave vectors n(a;) and #(x) and the combined system T{x) undergo changes determined by 
equations f^^^ - iP^^ . f^^^ - ipHi^ ^ and l|^^7H)l - l(2T7| . Besides the deformations induced by infinitesi- 
mal changes of the endpoints {aj}, all these deformations have the same general form, depending only on 
the function v{x) — 5V{x) that gives the infinitesimal deformation of the potential. We deal with them 
all on the same footing in the following proposition, which expresses the explicit form they take on the 
window r„(a;). 

Proposition 3.2 The deformation equations \2-7S\ ), \2-7S^ \2- 7?| ) are equivalent to the infinite sequence 
o/ 2 X 2 equations 

), (3-16) 

where the folded matrix of the deformation is defined by 
v{x)-\v{Q)n- 







-l,n-l 



+ In 



Ti— 1 , n 

VQv{x)n,n-l ^ Qv{x)nn 



(3-17) 



For the deformations in 
to changes in the potential 



id ^2-65\) - f2-69p . this gives the following equations corresponding 



hdcT^{x)=Cr;nix)Tn{x) (3-18) 

hdt^^,rn{x)=%.J,n{x)r„{x), (3-19) 

where the sequence o/ 2 x 2 matrices Cr-.n o,nd Tr,.j-n{x) are rational in x, with poles at the points {cr}, 
obtained by making the following substitutions in eq. ( {3-17}) . 

^ - tr,,j{x - Cr)' 





v{x) 


Tr,J ■ 


v{x) 


To,j : 


v{x) 


Tr,0 ■ 


v{x) 



-J-1 



,7=0 

— ln(a; — Cr) 



(3-20) 
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Proof. Using the definition (|2-7()(l of Vqu(x) and the extended ChristofFel-Darboux relation (|3-8|l . we 
have 



7nVQ«(a;)S„r = 7„ / dye-T^^fe) (v(^y) _ y(x)) n(y) 



n^(y)s^r(x) 

y-x 



+ (o, -e*^(^) I (^^^^^) n(y) 

= z;(Q)n°_ir(a;) - «(a;)n^ir(a;) - (0, VQi;(a;)i7(a;)) 
Applying the projector and noting that 

ilv{Q)Ul_,r = tlv{Q)-oT- 

we obtain 



1 / v{Q)n-l,n-l 

2 I 







(3-21) 

(3-22) 

(3-23) 
(3-24) 

(3-25) 



M,r„(x) = C^^.r(x) = C (^.r(:E) + (O, VQwC/(a;)) 
= tlv{Q)-oT{x) + {0,^lWQv{x)U{x)) 

— iu(Q)„_i^„-i 



z^i;(g)n°_ir(a:) - 

lni'^^^Qv{x)ilaTn{x) 



V{x) - ^v{Q)n 



9 '-'V^ yn— l,n— 1 

i 



(3-26) 
(3-27) 

r„ + (0,i^VQw(a;)J7(a;)) (3-28) 




w(Q)n,r: 



r„(x), 



(3-29) 



proving the relation (|3-17|) . 



Remark 3.1 Note that formula 13-1711 for the deformation of the measure in Proposition 13.21 as well as those 
below, 13-361 . 13-3711 . which are obtained through folding of the hd^ operator, could also be derived for arbitrary 
locally analytic potentials V{x), provided all the integrals involved are convergent |13l . However applicability 
of the subsequent isomonodromic analysis would be lost if the derivatives were not rational, since the resulting 
deformation equations would then have essential singularities. 

3.3 Folding of the endpoint deformations 

The case (12-2411 and (|2-65|) involving deformations of the locations of the "hard edge" end-points must 
be considered separately. 

Proposition 3.3 The following gives a closed system for the n-th window of eqs. {2-2^^ and \2-t)5]) 

h^aTn{x)=A^,n{x)Tn{x) (3-30) 

where 



ai — X 



V'n-l(ai) ipn-i{ai)ipn{ai) 



hKj 
' ~Y 
h>u_ 



-V-Liloi) 
V4(a.) 



(3-31) 



Proof. This is very similar to the proof of Prop. 13.21 Using the definition (|2-3()(l of the matrices Ai 
and the extended Christoffel-Darboux relation H3-8|l we have 



hdaTnix) = Hi'^daT{x) 



.i(y(^)-V(a.)) 



X — ai 



-n(a,) 
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flKi 



^i'l-M) 
V-^K) 



r„(x) 



rp^nfl-Xi , s T/ x„ ^ ^>^i 

ai — X I 



-^l-M) 







(3-32) 



Recalling the definition (|3-5|l of S„ and computing the matrix product yields the result in the statement. 
Q.E.D. 

3.4 Folded version of the recursion relations and hd^ relations 

We now consider the recursion relations H2-20I) . H2-63() and H2-75|l and the action of the Hdx operator in 
(|2-23|l . (|2-64|) and (|2-75|) which, in their folded form are given by the following. 



Proposition 3.4 The folded forms of the relations ^2-75}) and \2- 76] ) 

r„+i(a::) = i?„(x)r„(a;) , n > 1, 



where 



and 



with 



R„ 



X — Oi 



1 

7„+i 7„+i 



i/'„-i(ai)-(/'n(ai 



-V'Li(a») 
-■4'n-i{ai)'^n{ai) 



V'ix) 



V'{x) 





In 



(VQr(x))_,.„_, (VQr(x))„_, 

(VqV^'(^))„,„_i (VQr(x))„„ 

"(VQr(:.))„_,.„ -(VQr(x))„_,„_, 

(Vq^'(^))„„ -(VqV^'(:i:))„.„_i 





■ 


-7ri 




_7n 






Remark 3.2 Note that formula implies that 

Tt{V4x)) = V'(a;) . 

Proof: The folded form H3-33|) of the recursion relations follows directly from eqs. (I2-2U|) and 

XTTn{x) = 7„+i7r„+i(a;) + /3„7r„(a;) + 7„7r„_i(a;) 

x(t)nix) = 7„+i(/'„+i(a;) + Mn{x) -l-7„0„-i(x) + SnQ\/%eT^^^''^ . 

To prove H3-34|l . note that the folding relations H3-lt)|) may be expressed 
for any infinitesimal variation v = SV in the potential. Choosing 



K 



do 



,j+idtQ j '^to,idtoo, 



(3-33) 
(3-34) 

(3-35) 
(3-36) 



(3-37) 

(3-38) 
) 

(3-39) 
(3-40) 

(3-41) 
(3-42) 



J=i 
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we have 

Using H2-29|) and H2-64|) . we have 



V'{x) = 6V. 



pr+ , \7qV'{x)u ~hJ2 



^l.(V(x)-V(a,)) ^ 

Xi n(ai) 

X - a,. 



' ( gi(y(x)-y(a,)) ^ 

(v^'(Q)-o-5I^0r+ 0, Vq{5v){x)u- hJ2^, — — n(a,) 



4=1 



f ^j-{Vix)-V{a,)) 

(('^■^')(Q)-o-5I^*)^+ VQ(5F)(a;)J7- ?i^>^, — — n(a,) 



X — Oi 



r, 



(3-43) 

(3-44) 
(3-45) 
(3-46) 
(3-47) 



where we have used the deformation equations H2-24|l - H2-26|l . H2-65|l - H2-69ll . Applying the folded relations 
(EU, l|5T7|) and this gives 



where 



ai — X 



K 



K 



il^litti) -■il;n-i{ai)ipn{ai) 

-^n-lK) 

^2 (a,) 



(3-48) 

(3-49) 
(3-50) 



and 



Vnix) 



+ 



V'{x)^W'iQ)n-l..-l 



VQy(x)„_i,„_i VqF'(x)„_i,„ 
VqF'(x)„,„_i VQr(x)„„ 



"7n 
7n 



(3-51) 



It follows from H2-39|) that the diagonal V'{Q) terms in V„(a;) are cancelled by the sum in the last term 
of H3-48|l . giving the stated result 

Combining the differential, recursion and deformations relations (|3-34|) . (|3-33|) . (|3-18|) . H3-19|) and 
(|3-3U|I , the fact that the invertible matrices r„ are simultaneous fundamental systems for all these equa- 
tion implies the compatibility of the cross-derivatives; i.e., the corresponding set of zero-curvature equa- 
tions. 

Corollary 3.1 For n>0 the set of PDE's and recursion equations 

hdxTn{x) = I?„(a;)r„(x), hdaTn{x) = Ai-n{x)Tn{x) 

hdcTri{x) = Cr:n{x)Tn{x), Hdt^ j'^nix) = Tr^j.riix)Tn{x) 

T„+i{x) = R„{x)r„{x) ' (3-52) 

are simultaneously satisfied by the invertible matrices r„(x), and hence the zero- curvature equations 

[hdx - Vn, hdai - A;n] = 0, [hdr^ - I?„, Hda, - Cr-n] = 0, 
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[hda, - Ai-n, flda, - Tr.J-n] 
fidaiRn 



, [flda, = Ai-n, fidai " ^r-n] = 0, 



(3-53) 



are satisfied. 



Remark 3.3 (The Riemann-Hilbert method.) 

The Riemann-Hilbert method for characterizing orthogonal polynomials |10l provides an alternative ap- 
proach to deriving the results of this section. This is a well-established approach, and will not be developed 
in detail here, except to indicate briefly how it could be applied to deducing the differential and deformation 
equations satisfied by the f undamental systems. 

The fundamental system r„(a;) has, by construction, a jump-discontinuity across any of the contours defining 
the orthogonality measure, denoting the limiting values when approaching any of these contours from the left or 
the right byr)n,± we have the jump disconinuity conditions 







I 



(3-54) 



Furthermore, the local asymptotic behavior near the singularities at oo are specified as in section 1^?^ To be more 
precise the function Tn(x) has local formal asymptotic form, within any of the Stokes sectors. 



*(a;) 



Cr \ l + 0{x~Cr) e~2 



Cr 



Aj\l + 0{x~a,) e 



Co ( 1 + O ( - 

X 



-Xj ln(x — aj)(j^ 



^To(a;)<T3 + (n-E,, to,,.)o-3 ln(a;) 



(3-55) 



It follows from the usual argument based on Liouville's theorem that any two fundamental solutions two 
(with the same Stokes matrices, given in fact by the same matrices (13-5411 ) satisfying the above Riemann-Hilbert 
conditions are equal, within a constant scaler multiple. Also, from Liouville's theorem it follows that the first 
column of ^'(a;) consists of polynomials (the orthogonal polynomials). Using similar arguments one can show that 
the following matrix is rational with poles, of the correct order, at the singular points c,.,aj,oo: 



V,,(x) ■- a.*(s)*(s)'' + iv'(a;)l. 



(3-56) 



By comparing the local singular behavior of the logarithmic (matrix) derivatives of any two solutions and applying 
Liouville's theorem, it again that these globally combine to define rational matrix functions which give the 
deformation matrices with respect to the various parameters at the poles. 



4 Spectral curve and spectral invariants 

The aim of this section is to express the spectral curve of the ODE 13-341) (i.e., the characteristic equation 
of I?„(a;)) in terms of the partition function. In fact we will prove an exact finite n analog (Thm. I4.2|l 
of the formulae that are obtained by variational methods in the n ^ oo limit [S]. We start by expressing 
the explicit relation between the partition function and the spectral curve of the isomonodromic system. 
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4.1 Virasoro generators and the spectral curve 

To express the result in a compact form, introduce the foUowing local Virasoro generators 



VL''^ := J2 Mtr, 



d 



M=l 
do -J 



V 



(0) 

-J 



M+J 



M+.J 



dtr^M 

d 



dto. 



J^O,...,dr-l 



J = 0, . . . ,cfo - 1, 



M 



(4-1) 



(4-2) 



in terms of which wc define the following differential operator with coefficients that are rational functions 
of X 

_ V V - v^""^ ^ 



L ^ do — 3 

H.) :=V^^- V.T'^-(°^ 



— Cr dCr 

7_n . 1 T_n \ / r— 1 



(4-3) 



Theorem 4.1 The characteristic polynomial of the matrix V^ix) in the differential system is 
given by 



det (y - Vnix)) ^y^~ yV'{x) + ft^Tr 
= y^ - yV'{x) + hTin 



V'{M) - V'{x) 
M - X 

V'iQ) ~ V'{x) 
Q - X 



ln{Z„ 

^ X - a., 



1=1 



da, ln(Z„) 

T — n. 



d„-i 



= y^ - yV"{x) + n ^ toj+ix'^-^ - ^^©(x) lnZ„ 
and the quadratic trace invariant is 



J=i 



d„-i 



TiVnixf = V'{x)^ - 2n ^ to.j+ix^~^ + 2h^B{x)lnZ„ 



J=i 



Proof: The equivalence of H4-4f) and 14-5|) follows from the well-known relation 

(Tr(/(A/))) - Tr„(/(Q)) 



(4-4) 
(4-5) 
(4-6) 

(4-7) 
(4-8) 



for any scalar function f{x) for which the (Tr(/(Af))) is a convergent integral. The equivalence of (|4-6|l 
and (jlTjl follows from 

To prove 14-5|) we use the recursion relation H3-33|) and the explicit expression of I?„, we obtain 



Vn+l — RriVnRn ^ + HR' Rn ^ 



n' n -1 








7n + l 









Therefore, 



Tr (l?„+i(a;)2) = Tr (P„(a;)2) + 2hTT {Vn{x)Rn-^K) + h^Ti {R'^R 



= Tr {Vn{xf) - 2h 
= Tr (Vnix)^) - 2h 



V'iQ) - v'{x) 
Q - X 

V'iQ) - V'jx) 
Q - X 



}<^■4'lia^) 



(4-9) 
(4-10) 

(4-11) 
(4-12) 



L 



+ 2h^ V -^da, ln(/i„), (4-13) 
■^-^ X — a. 



i=l 



20 



where we have used eqs. (13-36(1 . ((3-37(1 in (14-12(1 and ((2-47(1 in ((4-13(1 . These equations imply that 



Q - X 



-1 L 



]3 



j=l i=l 



From the definition of Vi , we have 







00 ■ 




TTO 


4>0 






0i . 




. ^1 


01 



Using 1(3-12(1 . this gives 

det(2?i(x)) = h^-ile'^^^^^'^^l[ 

^ hohl 



ho 
ho 



V'{x) 



\v'ix) 



x ~ z 
X — z 



-dz 



-dz 



X ~ z 
-dz 



-dz 



1 



(r(:r)-r(z))V§(z) 



X — z 



dz + h^Y^ 



dz \x — z 

Ki^Po{ai) 
X — ai 



dz 



and hence 



Tt{VI{x)) = -2det(X>i(a;)) -f Tr(X>i(x))2 
^{{V'{x)f~2h 



{V'{X)-V'{Z))^1{Z)^ J£.^^^2(„^) 



= {V'{x) f - 2h 
Combining this with 1(4-14(1 gives 



X — z 

V'{x)^V'{Q) 
X - Q 



-dz -2h^Y 



Tr(I?„(x)2) = iV'ixJf - 2hJ2 ^^'^^^ " ^'^^^ 



i=i 



Q - X 



n L 



^-^ ^-^ X — ai 



33 



j=l i=l 



(4-14) 

(4-15) 

(4-16) 
(4-17) 

(4-18) 

(4-19) 

(4-20) 

(4-21) 
(4-22) 

(4-23) 
(4-24) 



which, taking the expression ((2-10() for the partition function into account, completes the proof of eq. ll4-5() . 

We now proceed to the proof of eq. ((4-6(1 . By expanding the third term on the right of ((4-5(1 . we 
obtain 



V'{x) - V'iQ) 



do-l 



J-1 



K d^ + l 



J-1 



^ J=l A/=0 

do-l do-l ,7-2 

^ <o,J+ia;^"'l+ E ^*'*o,J+i' 



1 



K dr + 1 



(x-c^)^-A^ 
J-1 



(Q - Cr) 



-A/-1 



1 



J=l 

do-l 



J=2 J\/=0 
do — 3 



r=l J=l 



M=0 



(x - c,.) 

d.. + l 



J-M W "^"^ 



-M-1 



J=l M=0 J=M+2 



iJ-M 



— ' ^—^ (x — Cr) 

= 1 J\/=l ^ ^ 



M 



tr,J-l{Q - Cr) 



M-J-1 



J=M 
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(4-25) 

Now recall that for any deformation matrix X corresponding to an infinitesimal variation d we have 

n—l , 

Y.X,,=--dHZ,,). (4-26) 
i=o 

Summing the diagonal terms of (|4-25|l up to 71 — 1 and substituting into (|4-5|l we therefore obtain 

do-l 

Tr(2?„(a;)2) = {V {x)f - 277;7 ^ h,j+ix-^~^ 

,7=1 

do— 2 do ,Q 1 /' '7 A 



M=l ,7=M+1 dto,J-M-l 

(91n(Z„) 



E E 7— Tw E (-^ - + : 



(x — Cr)^^ '' dtr 

2h^ J2 + Xi E (4-27) 

r—l j—0 2—1 



91n(Z„ 



Using eq. H4-26|l we finally get 

rfn — 1 rfo — 3 do — 1 

det(y-I?„(x)) = y2_y-^^/(^)+„;j^ ^^^^^^^J-i_^;j2 ^ ^ ( J - A/ - l)to,j+i 

J=l M=0 J=J\7+2 

+ E E E - + ^)^r,-.^p^ 

I 1 ain(Z„) f 91n(Z„) ^ ^^ ^^^ 

r—l i—l 

which completes the proof of eq. H4-6|l . 
4.2 Spectral residue formulae 

Theorem 14.11 which determines all the coefficients of the spectral curve as logarithmic derivatives of the 
partition function, may be expressed in another form, in which the individual deformation parameters, as 
well as the logarithmic derivatives with respect to them, may be directly expressed as spectral invariants. 
The characteristic equation of I?„ (x) 

det(y(a;)I-P„(x)) =0, (4-29) 

defines a hyperelliptic curve C„ as a 2-sheeted branched cover of the Riemann sphere, on which y is a 
meromorphic function. It follows from and Theorem l4.1l that y, viewed as a double valued function 

of X, has the same pole structure and degree as Vnix) at the points {cq — 00, Cr, a,}, but that the points 
{ai} are branch points. 

Let Y±{x) denoted the two values of y{x). Defining 
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it follows from the explicit expression (|4-5(l for the spectral curve that, near any of the poles co 
oo, ci, . . . ,ck, the two branches have the asymptotic form 



Y±ix) = ^V'ix)±^^iV'ix)y-W 



i \ , 1 / \ r ©(x-^-io-i) X 



(4-31) 



Theorem 4.2 The following residue formulte express the deformation parameters and the logarithmic 
derivatives of Zn as spectral invariants of the matrix I?„ (a;) . 

1 /■ Y~(x) 

to J = 77- f -^^dx, J = 1 ... do (4-32) 

'^^'^ Joe 

trj Cryy+{x)dx, r = l,...,K, J^l,...,dr (4-33) 

H^dt, In Z„=-^ <f r_ {x)dx = -nh, (4-34) 

h^dt, lnZ„ = --^ <f Y^{x)x''dx, J = 1, . . . , do (4-35) 

h^dt^,\nZ,,^--^ (f> Y^{x)- ^—jdx, r=l,...,K, J ^l,...,dr (4-36) 

ZlJTT J [X — Cr)'^ 

h^dc^ lnZ„ = -r^ / Y^{x)T;{x)dx, r = l,...,K, (4-37) 
S'^a^ lnZ„ = ^^£ Tr(I?2 (x))da: . (4-38) 

Proof. The second equality in eq. 14-34|l follows from the fact that io.o appears in the integral H2-1U|I 

defining Z„ only in the overall normalization factor e ~ . The proof of the other relations is based on 
formula (|4-27|l , which was proved in the demonstration of Theorem 14.11 From this formula it follows 
that, for all deformations except the ones with respect to io.doi ^o,do-i' ^-i' 

h'dlnZ,, = i res ^^Tr(I?2 (^)), (4-39) 

where dV{x) denotes the deformation d of the potential V, and the residues are taken at the corresponding 
singularity Cr- Considering the various deformations associated to the poles and end-points we have: 
At infinity 

^X..lnZ„ = l^i;e^|^Tr(P^(x)) 

^2/iin^r -2^Tr„^-^ + 2)^— -9.^.1nZ„^ 

- res #^Tr„m^ . ^Tr„g.^ , . ^ 1, . . . do - 2. (4-40) 
2=00 V (X) X — Q J 

Note that this computation does not provide the derivatives with respect to the two highest coefficients 
to^do and to,do-ij which will be computed below. Moreover we should remark that the last equality follows 
from the following interchange of order of integrals 



res ^Tr„m^™ = E - 44 / ^'^^^ ' ^'^^^ ^/(.)e-^-'-)d. 
x-Q ~^^=^V'{x)J^ x-z 



x=oo V'ix) „ ^ 
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The exchange is justified by the usual arguments observing that the expression ^ has no singu- 

larities at coinciding points x — z (away from the singularities of V). 
At the poles 



1 



h dt^ J lnZ„ = - res 



Tr(I?^(a;)) 



h^dr lnZ„ 



X^Cr V'[X) X ~ Q J 



(4-42) 



The last equalities in 14-42|1 are obtained by a similar argument used for the deformations at cq = oo 
here above. 

At the endpoints aj 



h^da^lnZn = i res Tr{Vl{x)) 

Z x—aj 

= i res f {V'{x)f - 2hTT, 

2 x—aj \ 

= ffda^ lnZ„ . 



V'{x)-V'iQ) 



2E' 



-da. lnZ„ 



(4-43) 



The relations 14-33|1 now follow immediately from (|4-31|l . The determination y_ has the asymptotic 
behavior 



1 



V'ix) 



x-Q 



O(x"2rfo-i) x^oo 
0((a; - c,)2'i'-+5) x^cr 

(4-44) 

Identities (|4-34|l . H4-36(l and 14-37(1 follow immediately from the expressions in ((4-4011 . ((4-4211 and the 
asymptotic forms (14-44(1 . as do the identities ((4-35(1 for J < do — 2. For the remaining two values of J 
(do — l,do)j we compute 

- res T^{x)Y^{x)dx = - res (y' {x) + 0{l/ x))—— 

x—oo x—oo V \X) \ 



w 



+ C'(x-3))j dx 

= - ^ fiTr„T;(Q) + ^"^a, lnZ„=?iTr„r^(g), (4-45) 



K 



where the last equality follows from eq. 12-391 (translational invariance). This identity, together with 
eqs. ((4-35(1 for j < do ~ 2 implies 



- res x'^°~^Y^{x)dx = 

x—oo 

which is the case J = do — 1 of ((4-35(1 . Similarly 



d 



dto,dQ- 



(4-46) 



K 



- res xmx)Y_ix)dx = - res {xV'{x) + V o + 0{l/x))-^ 

x=cyo x—oo ^ — ^ V (x) 



r=l 



w - 



+ 0{x-^))^ dx 
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K K L 

= n2fi2 _ nhY^trfi - ^ HTinQZiQ) + J2 h^ajda, lnZ„ 

r— 1 r— 1 J=l 

= fiTr„QT^(g) (4-47) 

where the last equahty holds because of dilation invariance. This, together with the above proves H4-35|l 
for J = d-Q. 

The last formula H4-38|) follows from equating the residues at the poles x = Uj in eq. 14- 7|) . 

Remark 4.1 In the formula l)4-;-{4ll . l|4-:-{5|l . l|4-:-i7|l . H4-H8II we may replace Y- by — 1/'(a;)); this corresponds 

to the fact that, in the large n limit, the behavior of Y{x) on the physical sheet (i.e. Y+) is related to the resolvent 
of the model by 

Y+ ^ V' (x) + h {Tr{x - M)'^) . (4-48) 



5 Isomonodromic Tau function 

5.1 Isomonodromic deformations and residue formula 

In this section we briefly recall the definition of the isomonodromic tau-function given in and compute 
its logarithmic derivatives in the present case in order to compare it with the partition function. This 
will lead to the main result of this section, Theorem l5.ll which explicitly gives this relation. 
Consider a rational covariant derivative operator on a rank p vector bundle over CP^ 

= a. - A{x) , (5-1) 

where the connection component A{x) is a p x p matrix, rational in x. Deformations of such an operator 
that preserve its (generalized) monodromy (i.e. including the Stokes' data) are determined infinitesimally 
by requiring compatibility of the equations 

d^'i'ix) = A'i'ix) (5-2) 
du.'i'ix) =U,{x)^ix) , i = l,... . (5-3) 

where in the second set of equations Ui{x) are also px p matrices, rational in x, viewed as components of 
a connection over the extended space consisting to the product of CP^ with the space of deformation pa- 
rameters {wi, . . .}. The invariance of the generalized monodromy of Vx follows 11 from the compatibility 
of this overdetermined system, which is equivalent to the zero-curvature equations 

[dx - A{x),du, - U,{x)] = , [du, - U^{x), du, - U,{x)] = (5-4) 

Near a pole x = of A{x) a fundamental solution can be found that has the formal asymptotic behavior, 
in a suitable sector: 

^{x) ^ C,Y,{x)e^'^'^^^ (5-5) 

where is a constant matrix, 

Y^{x) = l+0{x-c) (5-6) 

is a formal power series in the local parameter (x — Cy) (or 1/x for the pole at infinity) and Ty{x) 
is a Laurent-polynomial matrix in the local parameter, plus a possible logarithmic term tQ\ii{x — c). 
In the generic case Ty{x) is a diagonal matrix, and, more generally, may be an element of a maximal 
Abelian subalgebra containing an element with no multiple eigenvalues. The locations of the pols Cy 
and the coefficients of the nonlogarithmic part of Ty{x) are the independent deformation parameters. 
The deformation of the connection matrix A{x) is determined by the requirement that the (generalized) 
monodromy data be independent of all these isomondromic deformation parameters. 
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Given a solution of such an isomonodromic deformation problem, one is led to consider the associated 
isomonodromic r-function determined by integrating the following closed differential on the space 
of deformation parameters 

:= y res Tr iY-^Yl ■ dT^{x)) = dlnr-^*^ , (5-7) 

V 

where the sum is over all poles of A(x) (including possibly one at x = oo), and the differential is over 
all the independent isomonodromic deformation parameters , In the present situation A{x) is our 2x2 
matrix TDnix) and the (generalized) monodromy of the operator dx — 2?„(x) is invariant under changes 
in the parameters Cr, tr^j , aj and n. 



5.2 Traceless gauge 



For convenience in the computations we perform a scalar gauge transformation of the ODE by choosing 
quasipolynomials rather than polynomials. Explicitly we set 



^r.{x) := e-A^Wr„(a;) = 

An{x)=V,,ix)^^V'ix)l, 



iJm{x) 1pn{x) 



(5-8) 



where ^ 

V-n :=e-A^(^)Vn (5-9) 
In this gauge the matrix of the ODE is traceless and the infinitesimal deformation matrices are trans- 
formed correspondingly by addition of the identity element multiplied by the derivatives of — 57^^(2;) 
with respect to the parameters {cr,trj,aj}. This choice gives a consistent reduction of the general 

gl(p, C) isomonodromic deformation problem to sl(p, C). (To be precise, this would require a further, 

-1 1 

^-independent diagonal gauge transformation of the form diag(/i„^_j^, hn) to render the infinitesimal de- 
formation matrices also traceless.) 

At each of the poles cq :— 00, ci, . . . we then have the following asymptotic expansions. (To simplfy 
notation, the index n is omitted in labelling the fundamental system 5* and its local asymptotic form .) 



■^{x) ^ CrYr{x) exp 

Here we have set 

Yo{x) 1 > ; , Co 



CT3 



(5-10) 



1 + 

fe=l 




1 



\/hn-l 





7r„_i(cr)e 2ft [Cr - Q)n-i,oVhoe 2ft , 
7r„(c,.)e 2R (Cr - QjnoVhoe^^ 



(5-11) 



where Vr{x) = V{x) — Tr{x) is the holomorphic part of the potential at Cr- The asymptotic forms given 
by I|5-1Q|I - H5-11|1 follow from the fact that, in any Stokes' sector near c^, the second-kind solutions behave 
like 



'ipnix) 
^n{x) 



k=n+l 



^„(z)c-*^z'=-i = e^^(")x-"-i (|l + O (5-12) 



A-(^) / d/'^'^^^^-"(-)(l 



+ Oix-Cr)) 
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Near the endpoints aj we have 



(5-13) 



~ Aj ■ Yj{x) ■ exp [^>Cj \n{x — aj)(T_|_] , (7+ := ^ 
7r„_i(a,)e 2. e j^dz ^ 



1 




since the matrix 













X- z 



^„(aj)e e j^dz ^-^^ 



7r„_i(a;)e 2r e 2r J^dz ^IV — ^"^^ 



(5-14) 



7r„(a;)e 



(5-15) 



is analytic in a neighborhood of Oj and has the hmiting value indicated in H5-14(l . The function — J dz " 
in the exponential of the second matrix in this formula has the same singularity as >fj ln(a: — a^). (The 
signs in H5-14I follows from the orientation of the contour originating at flj).) 
The differential 15-7|) can now be written 



Mlnrf^ = -\Y res dT,(:r)Tr (y-iF/ag) + V res ^^Tr (Y^-^Y,' a 



2 ^ — ' x—c 

r=0 



where the differential involves the isomonodromic parameters only 



K 



r=0 



dt 



dCr 



r,J 



dc 



d d 



da 



(5-16) 



(5-17) 



We now derive residue formulae for the deformation parameters and the logarithmic derivatives of the 
tau function for our rational 2x2 isomonodromic deformation problem. These essentially are the same as 
the formulae of Thm. 14.21 giving the latter quantities in terms of logarithmic derivatives of the partition 
function of the matrix model^. 

Consider the quadratic spectral invariant near any of the singularities: by virtue of the asymptotics 
()5-10ll5^T^ we have, near Cr and aj respectively (setting 5* := 2nS -I- X^^Li ^i-,o) 



Tr(^2(^)) ^ ?i2r[,^((^/^-i)2) ^ ^2rp^. ^i^y-iyiy^^ ^ ^^rj.^ ^ Y.-^Y^as ( r; 

2 



SrClS 



T' 



SroS 



(5-18) 
(5-19) 

Taking the principal part at each singularity and using Liouville's theorem (since Tr^^ is a priori a 
rational function) we find 



TTiA'ix)) = h'Tvii^^-'f) = h'Tv ({Yf^Y;)') + i^Tr {Yr^Y;a 



K 



i:t{a\x)) = Y1 



T' 



SroS 



T' 



SroS 



Tr(y-iy>3; 



r,+ 



^The case of an arbitrary rank rational, nonresonant isomonodromic deformation problem will be developed elsewhere 
151 . together with further properties that allows us to view these as nonautonomous Hamiltonian systems, in which the 
logarithmic derivatives of the r-function computed below, are interpreted as the Hamiltonians generating the deformation 
dynamics. 
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E 



XjTT{Y'{aj)a+) 



(5-20) 



where the subscripts r,+ mean the singular part at the pole x = Cr (including the constant for x — cq — oo). 
Consider now the spectral curve of the connection hdx — A{x) 



w 



-TvA^ix) , 



(5-21) 



one immediately finds 



w±{x) — ± 



SroS 



, SroS 
X 



Tr(yr^y>3)) 



KjTT{Y'{aj)a+) 



Near any of the poles one has the asymptotic behavior 

1 



(5-22) 



1^/ SroS 

2 



2x ' 2«t;(.)V^'"*'*' 



©((a; — CrY'"^^) near x = c^- 



near x = oo 



(5-23) 



'>r,Tr(l^'(a,)a+) 



-(1 + 0{x — Gj)) near a; = aj 



y ?i(a; — aj) 
This immediately implies the following identities 

nh H — tr.o = =F res w±dx 



— ^0,,/ = T res — jw±dx , J > 1 

2 X — QO X 

ntr,j = T res (a; — CrY w±dx 

2 ' X — Cr 



and 



9 



Sir, J 

d 



9aj 



IM 



± res (a; — Cr) •'w±dx 

X — Cj. 

± res T^{x)w±dx 

X — Cr 

res (ii;±)^(ia; . 



(5-24) 

(5-25) 
(5-26) 
(5-27) 



In order to compare with the formulae given in Thm. (|4.2(l we note that the eigenvalues w of ^(a;) and 
Y of T>n{x) are related as follows due to the change of gauge H5-8(l 



Y± = ^V'{x)+w± . 
Comparing eqs. 14-35|I - H4-38|I with equations (|5-25|I - H5-27|I we obtain 



(5-28) 



■IM 



— res x-^V 

2j a;— oo 
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Z \ 1 

^^"^ = "27 .=1 r = 1, . . . , i^, J - 1, . . . , d, 

ft2ae,ln4^ = -i res T;(x)y'(x)dx, r = l,...,if, 
h^da^ In Z„ = In r^*^ . (5-29) 

These relations define a closed differential 

dln(^-^^ =:dln(^„) , (5-30) 

where the quantity Tn is determined up to a multiplicative factor independent of the isomonodromic 
deformation parameters {c^, Oi}j>i; but which may depend on n. This may be explicitly integrated 
to give 

l-(§)=i E -sWT,(.), (5-31) 

where we have chosen to include the integration constant Inn! for reasons that will be explained in the 
remarks below. Note that Tn does not depend on the end-points parameters {aj}, only those entering 
in the potential V . Since the definition of the isomonodromic tau function r,(^^ allows normalizations 
depending arbitrarily on the monodromy data we obtain the following result. 

Theorem 5.1 t/p to multiplicative terms that are independent of the isomonodoromic deformation pa- 
rameters {trj, Cr, aj}, the partition function = Zn[{tr^ji c^, aj}| [><]) of the generalized random matrix 
model and the isomonodromic tau function t^^'^ for the associated ODE are related by 

Zn — Tn -^nj (5-32) 

where Tn is given in \5-til}) . 

We conclude with two remarks regarding the multiplicative factor Tn relating Zn and r^^^. 

Remark 5.1 (Change of gauge) 

The factor T„ can be eliminated by another choice of gauge. The original gauge is such that near any 
singularity Cr 

*(a;) = aK(x)e-^^'-'")"=' . (5-33) 

(The point x = oo would require slight modifications in the argument which do not, however, change the result.) 
The original tau function satisfies 

h^d,r) Inr"^ = --resTriY^^Yla-idTr) . (5-34) 

We make a scalar gauge transformation, depending on an arbitrary parameter c, of the form 

* := *e-t^("' = a?;e-(5"^ + ^^)^'- (5-35) 

where 

Yr ■- F,(a:)e-^^''-<"^ (5-36) 
and fr{x) is the regular part of V{x) at the point Cr 

fr{x) ■.= Y.T,{x) . (5-37) 
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The tau function for the gauge transformed system satisfies 



d^r) In rl"' = res Tr (y-Wr (a-i - 1 1) dT,) (5-38) 
= resTr [(y-^F,! - |/;) (^3 - |l) dT.] (5-39) 



= — res /;drr -I- res TTiYr'^Y^aadTr) res Tr(y^" V;) dTr (5-40) 

= — res frdTr + f^(r) Inr^^^ , (5-41) 

2 Cr 



where we have used that fact that 
RecaU that 



(5-42) 



d(r) In tI^^ = du-\ In Zn d(r) In J-n ~ di^-. In Zn H — res f^dT,- (5-43) 

2 2 



so that 



d(r) Inr""^ = — - — d(r) J"„ + d(r) In 2„ . (5-44) 



It follows that if c is chosen so that c +1 = (i.e. c = ±j) we have 

up to a multiplicative factor independent of the isomonodromic deformation parameters. 

Remark 5.2 (Schlesinger transformations) The shift n n + 1 given by the ladder matrix Rnix) corresponds 
to an elementary Sc/iZesmger- transformation at infinity. (For more details on Schlesinger transformation see I12|.'l 
In general a Schlesinger transformation corresponds to a shift by integers in the spectrum of the logarithmic 
terms of the matrix T[x) entering the formal asymptotics. We see from formula 15-101 1 that the shift n ^ n + 1 
amounts to increasing the first diagonal entry of the logarithmic term by 1 and decreasing the second entry by 1. 
As shown in |I2] the two corresponding tau functions are related by 



IM 

= (yo;i)i2 (5-45) 



Tn + l 



where the right hand side is the (1,2) matrix entry of the matrix Yo^.i in eq. 1)5-11^ . A simple computation (using 
eq. (I5-I2t ') shows that 

(yn,0;l)l2 = /in, (5-46) 

in agreement with Thm. 15.11 since (by eq. 12-1011 ) 

%ti = ;i„(n + l) . (5-47) 

This explains why the integration constant Inn! has been included in eq. 15-3111 : it assures that the dependence on 
the discrete parameter n in the relation between Zn and t^*^ is consistent with that implied by the Schlesinger 
transformations. 
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